Simulating and Detecting the Quantum Spin Hall Effect in Kagome Optical Lattice 
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We propose a model which includes a nearest-neighbor intrinsic spin-orbit coupling and a trimer- 
ized Hamiltonian in the Kagome lattice and promises to host the transition from the quantum spin 
Hall insulator to the normal insulator. In addition, we design an experimental scheme to simulate 
and detect this transition in the ultracold atom system. The lattice intrinsic spin-orbit coupling is 
generated via the laser-induced-gauge-field method. Furthermore, we establish the connection be- 
tween the spin Chern number and the spin-atomic density which enables us to detect the quantum 
spin Hall insulator directly by the standard density-profile technique used in the atomic systems. 
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I. INTRODUCTION 

Optical lattice system has gradually become a promis- 
ing platform to simulate and study a lot of quantum 
phenomena in condensed matter physics because almost 
all parameters of the system can be well-controlled [lj. 
The recent theoretical and experimental progress in laser- 
induced-gauge-field [2|-[7j makes it a hot spot to study 
topological quantum states [8[ in cold atoms system 1 9- 
fl5j . Subject to the compounds' natural properties [16l |, 
the famous Haldane model [17] proposed two decades ago 
has not been confirmed by experiments because the re- 
quired periodic magnetic field cannot be easily imple- 
mented in actual material. For topological insulator, an 
important extension from the Haldane model to a time- 
reversal invariant system |18l423J , only a few materials are 
confirmed currently to have such exotic topological prop- 
erties in the nature |24J27j , because the existence of such 
properties require relatively strong spin-orbit (SO) cou- 
pling. However, in optical lattice system, we can engineer 
the lattice Hamiltonian to guarantee that this system 
hosts these novel topological phases [l4|. Up to now, the 
neutral-cold-atom integer and fractional quantum Hall 
effects have been studied [28-30] and also the realization 
of Haldane model has been designed by using the laser- 
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induced-gauge-field method in optical lattice 

In this paper, we propose a scheme to simulate and 
detect the 2-dimensional (2D) quantum spin Hall (QSH) 
insulator in a Kagome optical lattice with a trimer and a 
nearest-neighbor SO coupling term. With laser-induced- 
gauge-field method, one can design a variety of lattice SO 
couplings [5|, I31M33J , which is convenient for us to study 
the 2D topological insulator in the optical lattice J9J-KL1J . 
However, the original proposal [18j of realizing 2D QSH 
insulator in honeycomb lattice requires the next-nearest- 
neighbor hopping amplitude. This requirement causes 
doing experiments in the optical lattice difficult because 
high barrier makes the next-nearest-neighbor tunneling 
very small. Recently, a model [34] raised in Kagome lat- 
tice also requires the next-nearest-neighbor hopping. In- 
terestingly, there is another QSH insulator model in the 



complicated Kagome lattice. The inspiration comes from 
the fact that spin chirality in ferromagnetic Kagome lat- 
tice exerts important effects on orbital magnetic moment 
and anomaly quantum Hall effect [35]. We find that a 
Kagome optical lattice with the trimer and SO coupling 
terms can host the 2D QSH insulator phase with only the 
nearest-neighbor hopping. Since the model only involves 
the nearest-neighbor hopping, it would be easier to be 
implemented in cold atomic experiments. 

Furthermore, compared to the condensed matter sys- 
tems, we find that detecting QSH insulator has more 
advantage in optical lattice system. Due to the time- 
reversal symmetry, Chern number can not be taken as 
a topological invariant to characterize the QSH insula- 
tor in real electron system [36|. Chern number Cf=+1 
(Ci=— 1) in a QSH phase for the up-spin (down-spin) 
electrons and the total Chern number C=Cf+C^=0. 
Since one can not distinguish the contributions of the 
conductance from the up or down spin electrons at cur- 
rent technology, it is impossible to determine whether the 
system is in the QSH phase or normal phase by measur- 
ing Hall conductance. However, in cold atomic systems, 
it is the atom's internal states that represent the spin, 
not real spin, which brings certain benefits to measure- 
ment. One can directly measure the spin Chern number 
to determine whether the system lies in the QSH phase 
or not because optically measuring the atomic internal 
states is very simple. In this paper, we demonstrate that 
the method developed to detect the Chern number in 
cold atomic systems [12|, [30| can be put forward further 
to measure the spin Chern number, and thus we establish 
the connection between the spin Chern number and the 
spin-atomic density which enables us to detect the topo- 
logical Chern numbers directly by the standard density- 
profile technique used in the atomic systems. 

The paper is organized as follows: In Sec. II, we intro- 
duce the model in Kagome lattice with both trimer and 
spin-orbital coupling. This model can realize the QSH 
phase with only the nearest-neighbor hoping terms; in 
Sec. Ill, we explain how to simulate this model in cold- 
atom optical lattice, which includes designing the lat- 




FIG. 1. (Color online) (a) Schematic picture of the nearest 
neighbor intrinsic SO coupling in 2D Kagome lattice. The 
up-spin atoms hop along (against) the arrowed direction with 
amplitude iXso (— iAso). For the down-spin atoms, the ar- 
rows are reversed. The dashed line represents the Wigner- 
Seitz unit cell, which contains three independent sites (A, B, 
C). (b) the trimer Kagome lattice. Hopping amplitude corre- 
sponds to t+K (t—K) for the thick (thin) bonds, (c) Contours 
of the effective magnetic field for up-spin atoms defined by 
Eq. ([21]). 



tice SO coupling by using the laser-induced-gauge-field 
method; we present the method to detect the QSH phase 
in Sec. IV and give a brief summary in Sec. V. 



II. MODEL 

Let us consider the tight-binding model for two- 
component fermionic atoms on the Kagome optical lat- 
tice, which consists of three triangular sublattices A, B 
and C (Fig. 1). The spin-independent part of the Hamil- 
tonian is given by 

-tio = to / J ^m,n,a a rn,n,a~^brn-l,n,a a rn,n,a 
mnoi 

i,n,a C m,n,a+am jn +i ia Cm,n,a J +H.C, (1) 



-ha, 



where to is the hopping amplitude between the nearest 
neighbor link, (m, n) labels the Kagome unit cells with 
the unit vectors bi = (2,0)a and b2=(l, VS) a, a^ n a 
(a>m,n,a) is the creation (annihilation) operator of an 
atom with spin a (up or down) on lattice site (m,n) on 
sublattice A (an equivalent definition is used for sub- 
lattice B and C). For simplicity, we choose to=l as the 
energy unit and the distance between the nearest sites a 
as the length unit throughout this paper. 

By using the Fourier transform of atomic operators 



am,n ' a ~ Vn ^ 



— ?'k U A 



(2) 



the Hamiltonian (JTJ can be diagonalized in the momen- 
tum space as 



ffo = XX(fto(k)®I 2 x2)tfk, 



(3) 



where I2 X 2 is a 2x2 unit matrix. The six-component 
basis operator ^k=Okt, &kt > c kt , a k ;, 6 k +, c k ;) T describes 
an atom in the corners of the Kagome lattice unit cell 
(Fig. 1) with spin a (*f and \). Here Ho(k) is a 3x3 
spinless matrix given by 



2 cos Pi 2cosP 3 
H (k)= [ 2 cos Pi 2cosP 2 
2cosP 3 2cosP 2 



(4) 



where we have defined Pi=k-ai, P2=k-a2, P3=k-a3 with 
ai=(l,0), a 2 =(-l/2, v / 3/2), and a 3 =(-l/2, -V5/2) 
representing the displacements in a unit cell from A to B 
site, from B to C site, and from C to A site, respectively. 
In this notation, the first Brillouin zone is a hexagon with 
the corners of K=± (2tt/3) a x , ± (2tt/3) a 2 , ± (2tt/3) a 3 , 
and two of which are independent. 

The energy spectrum for spinless Hamiltonian Ho(k) is 

characterized by one dispersionless flat band (e lk =— 2), 
which reflects the fact that the 2D Kagome lat- 
tice is a line graph of the honeycomb structure 
[37J, and two dispersive bands, e^L^=l^f\^4b^ — 3 

with bk=^2 i=1 cos 2 (k • a^). These two dispersive 
bands touch at Dirac points K and exhibit a cusp, 
6$ )k =(lTV3|k-K|). 

When the following intrinsic SO coupling term is taken 
into account in the Kagome lattice model, as illustrated 
in Fig. 1(a), the gap will be opened at the two inequiv- 
alent Dirac points. The tight-binding expression for this 
SO coupling Hamiltonian can be given as follows: 

#SO=2Aso 2^ ( &m,n cr * a m,n+& m -i ) n cr * a m,n 

m,n 

+ C m,n a ^m,n+C m+1)n _ 1 (7 z m)n 



+a 



i,n cr 2 C m,n _ r tt 777, ?n +i cr 2 C m,n ) +-H..C, [b) 



where Aso is the SO coupling constant, <j z is the Pauli 
matrix and aJn,n = ( a L,n,t> a ln,n,i)' Taking the Fourier 
transform ([2|) and considering the ^ k below Eq. (|3|), we 
have #so=Z) k V ; k / ^so(k)^k, where 



Hso(k) 



with 



H±(k)=±2i\ so 




H+(k) 
H-(k) 



— COS Pi COS P3 

-cosfb 
cos P 2 



(6) 



• (7) 



This SO coupling destroys spin SU(2) symmetry and 
opens a band gap Aso = a/3 Aso at Dirac point. 

Lattice trimerization can break inversion symmetry of 
Kagome lattice and also open a gap at Dirac point |34j . 



It is described by 



-"trim— / J 1^ ( bm,n,a a m,n,a ^m-l,n,a a m,n,a J 

mna 

+ ft \C-rn ,71, a,brn, n , a ~ C m+l,n-l,a' :) rn,n,a J 



x ra,n+l,a ?™,n,a 



H.c, (8) 



where ft describe an alternating pattern of bond hopping 
integrals along the three principal spatial directions as 
illustrated in Fig. 1(b). Taking the Fourier transform 
(|2|) again, the trimerized Hamiltonian can be rewritten 

as ^ t rim=Z)k^kC^trim(k)(g)l2x2)^k with 

—ft sin Pi ftsinP3 

H t rim(k)=2z( ft sin Pi -ftsinP 2 ) (9) 

-ftsinP3 ftsinP2 

for both spin components. 

We take above SO coupling and trimerized Hamilto- 
nian as perturbation, which means Aso^ an d ft<Ct. Al- 
though both of perturbations can bring gaps at Dirac 
points independently, these two gaps have different topo- 
logical nature. As we will show in the following, the 
former is non-trivial and quantum spin Hall effect will 
occur if Fermi energy level locates in the gap; the latter 
is a trivial gap. 

In order to prove above assertion, we expand the total 
Hamiltonian % (k)=Ho (k)+Hso (k)+%trim (k) at two 
inequivalent points K± = (±^,0), then take kp per- 
turbation theory to get its effective Hamiltonian. At last, 
projecting it onto bands 2 and 3 subspace, we get four 
independent Dirac Hamiltonian, 



n 



K . 



SVFk x T z + SVpkyTx 



(10) 



where s=±l and <j=±1 represent different valleys K± 
and spin indices, respectively, vp=y/^t is Fermi velocity, 
and Ti are Pauli matrices with i=x, y,z. The Dirac mass 
m sa =\^3a\so — 3sft. 

We now turn to address the topological properties of 
this model. The system is in the QSH phase when the 
SO coupling dominates in the condition |Aso| > |>/3«|; 
otherwise, the system is in the normal phase. It can 
be proved by directly calculating the Z2 topological in- 
variant [13, [20j . Alternatively, one can also understand 
the topological phases from the view of spin Chern num- 
ber |36j. In order to do that, we apply a uniform mag- 
netic field along z-direction with gauge vector A (r) = 
(0, Sx, 0) and obtain the Landau energy levels 



K a = { 



m sa sgn (eB) 



±^2nhv' 2 F \eB\ 



- mi 



n = 

n = 1,2,3, • • • 



(ii) 



By using the Green function theory [12[, we can get the 
Chern number 

C sa = ^sgn(£ s<T ) = i S gn(m sff sgn {eB)) (12) 
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FIG. 2. (Color online) Illustration of the effect on Chern 
number by the sign of n=0 Landau energy level. The left 
figure corresponds to disperse relation of Eq. (|1Q|) without 
magnetic field. The right two figure show the Chern numbers 
when magnetic field is applied and Fermi energy locate at 
different interval for two cases with Eq > and Eq < 0. 
Here Eq denotes the n=0 Landau energy level. The dashed 
line represents the zero energy and the doted lines give the 
original band gap. 



between the energy interval — \m sa \ </i< \m sa \ with \i 
being the Fermi energy (Fig. 2). Actually it is the 
sign of n=0 energy (zero-mode) which determines the 
Chern number. As an example, We focus on /i=0 
and take eB>0 throughout this paper to see the dif- 
ference more clearly between QSH phase and normal 
phase. For QSH phase, we take Aso>>/3^>0, and 
it is easy to see that, for up-spin atoms, ra s= ±i^>0 
and then C^=C+i j ^+(7_i j ^=1. For down-spin atom, 
m s =±i,i<Q and C;=C + i^+C_i^=-l. However, if we 
take v // 3^>Aso>0, which corresponds to the normal 
phase, C^=C^=0 for both up- and down-spin atoms. 
In the QSH phase, the Chern numbers with different 
spin components have same value but with opposite sign. 
Whereas if it is in normal phase, the Chern numbers 
equal to zero for both up- and down-spin atoms. 

To further understanding the topological properties of 
the model, we show the edge state effects in Fig. 3. From 
Fig. 3(a), we can see that there is a pair of chiral gapless 
edge states for every band gap when the SO coupling 
dominates. This means that the system is in topologi- 
cal insulator phases at 1/3- and 2/3-filling. When only 
trimer term exists, it opens a band gap at Dirac point but 
no edge states connect the upper and lower bands (See- 
ing Fig. 3(b)), therefore the system is in normal insulator 
phase at 2/3-filling. On the other hand we also see that 
the trimer term cannot open a gap between the band 1 
and 2 (Fig. 3(b)). Therefore the system at 1/3-filling 
will be still in the topological insulator phase when two 
perturbations are present but the trimer term dominates 
(Seeing Fig. 3(c)). 




FIG. 3. (Color online) The band structure of the lattice model 
in the stripe geometry. We take Aso=0.1, k=0 for (a), Aso=0, 
k=0.1 for (b) and A S o=0.05, k=0.1 for (c). 



III. SIMULATION 

In this section we introduce an approach to simulate 
the Kagome lattice with the trimer and SO coupling 
terms in an optical lattice system. To this end, two 
problems need to be solved. One is how to generate 
the Kagome optical lattice with the trimer terms; the 
other is how to simulate the lattice SO coupling. As for 
the first problem, we use the superlattice technique ad- 
dressed in Refs. [38l-[40| , that is, three super-laser beams 
with the same wave vector length but different polariza- 
tions are applied along three different directions: e y and 

^e^i^e^, respectively. In our proposal, each super- 
laser beam consists of four large detuned standing- wave 
lasers with the same polarization but different wave vec- 
tor length in the x-y plane. The total potential is thus 
given by 

3 

V (r)=Vb Y^ I cos Orr+3<W2) + 2 cos (k r r/3+<^/2) 



+4cos(ki-r/9+£i<p/6) 
+Ccos(ki.r/9+^ (y?/6+7r/2))] 2 



(13) 



with the wave vectors ki = (^, ^)fc, \a 2 =(— ^, \)k, 
k3 = (0, 1) k and 61=62=— 63=1. Firstly, we consider the 
case with (=0. One can get a triangular lattice when 
cp=0 or 2tt and a Kagome lattice when 0<cp<27T. A uni- 
form Kagome lattice corresponds to ip=7r, as shown in 
Fig. 4(a). When cp takes other values, one can obtain 
the trimerized Kagome lattice accompanying distortion 
of the lattice structure [38l-l40| . When we increase the 
strength of trimerized Hamiltonian, the Kagome lattice 
will be distorted. To overcome this defect, we add an- 
other laser beam which corresponds to the £^0 in the Eq. 
([T3]) and assume ip=7r all the time. The added laser will 
interfere with primary lasers and generate the trimerized 
Kagome lattice. ( is an adjustable parameter to con- 
trol the strength of trimerized Hamiltonian. With this 
method, the Kagome lattice will not have obvious offset 
from the uniform one even if £ takes a relative large value. 
As shown in Fig. 4(b), the parameter (=1-5 is chosen as 
a typical example. 

We now focus on how to simulate the lattice SO cou- 



pling. Using the laser-induced-gauge-field method, it was 
proposed that both Abelian and non-Abelian gauge fields 
can be simulated in cold atomic system. In addition, the 
experiments to achieve such artificial gauge fields have 
been reported [6|, [7J] • Interestingly, it was proposed that 
a periodic magnetic field, which is not easy to be real- 
ized in a condensed-matter system, can be created by two 
opposite-traveling standing- wave laser beams [5[. In the 
following, we will show in an explicit manner that the 
artificial gauge field proposed in Ref. [5] is equivalent 
to a SO coupling. Furthermore, the required lattice SO 
coupling addressed in the previous section can also be 
achieved in a suitable configuration of the laser beams. 

In optical lattice systems, cold atoms can hop between 
adjacent sites. According to Peierls theory, the additional 
gauge vector potential A makes the hopping obtain a 
phase factor exp (i ~ J A ■ dl) , where the integral along 
the hopping path. If the atom has multiple states, which 
correspond to different spin components, vector poten- 
tial A should be a matrix. After taking such an approx- 
imation, that is t?r=tij, here a, f3 on behalf of any spin 
index, the correction to atom hopping between different 
sites coming from gauge field is equivalent to a unitary 
operator [41[. For a two-component atom system, the 
unitary operator can be written as 



Ui, 



cosc^ 



i&ij sin c^ 



(14) 



where z, j represent different site indices, a^- is the gauge 
flux and depends on the hopping integral. Here aij is the 
Pauli matrix, whose specific form depends on the gauge 
vector potential. Therefore, a tight-binding Hamiltonian 
of atoms can be written as 



H 



(id) 



Zijdj U ij&i 



R.c.) = H + H so 



with 



H = ^2 (Uj c °s OLij) ajjCLi H- H.c. 



(15) 



(16) 



(iJ) 



(iJ) 

Here (i,j) denotes the nearest-neighbor hoping, ai (a\) 
the creation (annihilation) operator on site i. The first 
term is the normal Hamiltonian, while the second one is 
equivalent to SO coupling. From above equations, we can 
observe that: (i) By choosing suitable Pauli matrix <j^-, 
one can simulate various SO coupling existing in actual 
materials. For our model, we should choose <Jij=cr z . (ii) 
Through adjusting the gauge flux a^, one can change 
relative strength between the two terms. For example, 
otij=0 corresponds no SO coupling, while aij=ir/2 is 
equivalent to only the SO coupling interaction existed 
in the system. This facilitates us to study the nature 
brought by SO coupling. 

We now turn to address the approach to achieve the 
intrinsic SO coupling described in Eq.(|5j). To this end, we 
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FIG. 5. (Color online) (a) Illustration of the three-level A- 
type atoms coupling with two laser beams with Rabi frequen- 
cies Q p and Q c . (b) Another different hyperfine level \F) for 
detection which can be coupled to the excited state |3) by 
Rabi frequency Qf- 



FIG. 4. (Color online) (a) A uniform Kagome lattice for <p—ix 
and C— 0- (b) A trimer Kagome lattice for ip—ix and £=1.5. 



consider a cold atomic system with each atom having an 
A- type level configuration (Fig. 5(a)). The ground states 
|1) and 1 2) are coupled to the excited state |3) through 
spatially varying standing- wave laser field, with Rabi fre- 
quencies Q p =Qsm0e~ lSl and Q c =Qcos0e~ lS2 , respec- 
tively. With rotating-wave approximation, the laser- 
atom coupling Hamiltonian is given by 



H h 



Q p 

- I o o n c 
nt n* -2A 



(18) 



with the eigenstates (the dressing states) 



|X1> = 



iSi 



\D-e- 



iS 2 



sin 6 1 2) 



IX2) = cos ip sin 0e~ iSl 1 1) + cos cp cos 0e~ iS2 |2) - sin cp |3) 
|X3) = sin ip sin 0e~ iSl 1 1) + sin cp cos 0e~ iS2 |2) + cos ip |3) 



and eigenvalues A n =i 5 2,3=0,| (A^yA 2 +£l 2 ). Here, 
single-photon detuning A=ujs— ui— uu p , with CJ3, 
lji, u p the intrinsic frequency of atom states |3), 

In the new ba- 

IX2), |X3)}, the primary 

H=^+H int (r)+T> (r) can 

H=^j (-ihV-A) 2 +V with 
V being matrix 



|1) and laser fi p , respectively. 
sis space |x)={|Xi) 
atom Hamiltonian 

be rewritten as 
M the atom mass, 



_ p 

~2M 
~~2M 

A and 



with matrix element 



, n ^m=ih (Xn (r) I VXm 0)), 

V ni m=X n (r) S n ,m+(Xn (r) | V (r) \xm (r)), respectively. 
One can see that in the new basis the atom can be consid- 
ered as moving in gauge potential A, which corresponds 
to an effective magnetic field B e //=(V x A) — 4 (A x A) 

am. 

We focus on the subspace spanned by the two 
lower eigenstates {|xi) , IX2)}, which is redefined by 
Ixt) = Ixi) and \xi) = |X2> with the spin language. This 
gives an effective spin- 1/2 system. For the large detuning 
(A^>ft) case, both states |xt) an d Wi) are stable under 
atomic spontaneous emission from the initial excite state 
1 3). Furthermore, we assume the adiabatic condition, 
which requires the off-diagonal elements of the matrices 



A are smaller than the eigenenergy differences |A m — A n | 
(m, n=l,2,3) of the states \xm)- Under this adiabatic 
condition and in the basis space {|xt) > 1X4.)}? the gauge 
potential A becomes diagonal and takes the form [5j 



L t 







A. 



(19) 



with 



L t 



-A; = ft (VSi cos 2 0+VS 2 sin 2 0) 



Here we neglect the correction to nearest-neighbor tun- 
nelling brought by the change of potential V (r) because 
of the large detuning approximation. 

We consider a specific configuration of the 
laser beams with two opposite-travelling standing- 
wave laser beams [5|, LL2J? which take the 
Rabi frequencies 



(]p=Hsm 



-i(k±x+k z z) 



fay+j) e ^ klX+k ^ 



The 



and Q c =Qcos(k2y-\-j) 
effective gauge potential is generated as 
A^=— A|=ftsin(2/c22/+2</>i) fae x +k z e z ). Here 

ki=ksinOi cos 02, k 2 =k cos 0i , and k z =kshi6i sin 02 
with k the wave vector number of laser, 0\ the angle 
between the wave vector and e y axis, 02 the angle 
between e x axis and the plane consisting of wave vector 
and e y axis. We emphasize that the choice of wave 
vector &2 °f the laser beams must be a multiple of 
tt/^3 in order to be commensurate with the optical 
lattice. We take k2=7r/V3. The Peierls phase factors 
for the nearest neighbor hopping in Fig. 1(a) are 



.y/Skj. 



„fci 



with a=±l representing the up- and down-spin. Con- 
sidering the symmetry of Kagome lattice, the vector 
potential A is rotated by ±27r/3 to obtain the other two 
vector potentials. Therefore, the total effective vector 
potential and magnetic field can be written as 



Aeff=ahki 



( sin (2&22/) — cos fay) sin( V^/^x) ) 



Vr 



-y/3 sin fay) cos(V3k2x)e y 



-a — -=- [2 smfay) sin (irx) 

V3 



(20) 



■cos(2k 2 y)}e z . 

(21) 



It should be noticed that we have dropped the e z com- 
ponent in Eq. ([20]) because the integral for Peierls phase 
is only in x-y plane. The contours of the magnetic field 
for up-spin are plotted in Fig. 1(c). However, the to- 
tal accumulated phases for the nearest-neighbor hopping 
along the arrowed directions in Fig. 1(a) are 



^61=^45 = ^34 =( Pl2 = ^23=^56 = 



2 k sin 6 1 cos #2 

-a =aip. 



(22) 
We must retain 6\ to satisfy /c2=7r/v / 3. However, 
we can alter cp, which controls the relative strength 
between SO coupling and ordinary hopping terms, 
by changing the angle #2 in the x-z plane. So, 
we can replace Uj cos a^ —>t cos cp in Eq. ([T6]) and 
Uj sin aij—>t sin (p^Xso in Eq. ([T7j) and therefore get the 
intrinsic SO couping model in our cold-atomic Kagome 
optical lattice. 



IV. DETECTION 

For quantum Hall effect in two-dimensional electronic 

gas in condensed matter system, Hall conductivity a xy 

2 
and Chern number satisfy the relation <j xy = ^-C with h 

the Planck constant and e electronic charge. Thus one 
can usually detect the Chern number through measur- 
ing Hall conductivity. However, the detection of the spin 
Chern number in actual material system is challenge be- 
cause one can not distinguish the contributions of the 
conductivity from up- and down-spin electrons, respec- 
tively. 

Unlike in electronic system, we will show that a sig- 
nificant advantage of atomic system is that the spin 
Chern number can be directly verified by using the sim- 
ilar method to detect the (mass) Chern number. It has 
been shown that the conductivity a xy (Chern number) 
is related to the atomic density from the Streda for- 
mula <j xy =dp/dB I M5 t when a uniform magnetic field B 
is applied in the system. Thus one can measure the 
Chern number through the detection of the density pro- 
file, which is a standard detection method used in atomic 
system [12, [30j . Since the internal state-dependent im- 
age has also been achieved, the similar method can be 
straightforward expanded to measure the spin Chern 
number. Therefore we may establish the relation be- 
tween the spin Chern and atom density, enabling us to 
detect the OSH phase of the system via density-profile- 
measurement technique. 

We now introduce how to detect the QSH phase in 
ultracold-atomic optical lattice. In cold atomic system, 
the informations of quantum states are usually measured 
from the density profile image. So we will first derive 
spin-atomic density from the Dirac Hamiltonian ([TO]) , 
and then obtain the information of the Chern number 
of the system. 

The uniform magnetic field can be simulated by rotat- 
ing the optical lattice at a constant frequency u=eB/2M. 
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FIG. 6. (Color online) The spin-atomic density in units B/fo 
as a function of the normalized chemical potential p/mo with 
the definition hvp \eB\=4\mo\ . (a) and (b) correspond the 
QSH phase case with v / 3A 1 so>3^>0; (c) and (d) correspond 
the normal phase case with V3\so<3k<0. p^ (pi) denotes 
the up-(down-)spin-atomic density. 



From the Landau energy levels obtained at Eq. ([TT]) . we 
get the density [12] for every Dirac Hamiltonian ([TO]) as 



+00 

n=0 

B m s 

20o \m 3 



sgn (p) B 



20o 



[9(M-£ n+1 )+e(M-£ n )] 



r e(|ra sa |-|/i|). 



(23) 



where 9 stands for the unit step function, p is the 
chemical potential and 0o the flux quantum. The sec- 
ond term of Eq. ([23]) is the atom density coming from 
n=0 (zero- mode), similarly the discuss of the Cherm 
number (Seeing Fig. 2). The calculated spin-atomic 
density p a =pj ra J rp- a (a=t,i) in unit of B/<po is plot- 
ted as a function of the normalized chemical potential 
p/mo (for hvjp \eB\=4 \itiq\ ) in Fig. 6. It is essen- 
tial that the spatial density profile is uniquely deter- 
mined by the function p (p/mo) in the local density ap- 
proximation, which is typically well satisfied for trapped 
fermions. We focus in the point p=0 where the Dirac 
mass makes its central effect on the density or Chern 
number. When the system is in QSH phase, as an exam- 
ple, we take the parameters >/3Aso>3^>0, specially as- 
suming |m + ^| = |m_||=0.5 |m_^|=0.5 |ra+^|=0.5rao, then 
we have ra_^>ra+^>0>ra_^>ra+^. The spin- atomic 
density is shown in Fig. 6(a) and (b). It is easy to 
see that the up-spin-atomic density p^=B/(j)o>0 (Fig. 
6(a)) and the down-one p\=— 6/</>o<0 (Fig. 6(b)). 
For the normal phase, we take 3^>v // 3Aso>0 and as- 
sume |m+^| = |m_||=0.5 |m_^|=0.5 |ra+^|=0.5rao, which 
means ra_^>ra+^>0>ra_^>ra + ^. We obtain p^,i=0 
(Fig. 6(c) and (d)). The spin- atomic density shows 
the similar relation comparing to the conclusion of spin 
Chern number theroy. This is not surprising because ac- 



cording to the Streda formula and from Eq. (j23]h it is 
easy to obtain the relation between the spin Chern num- 
ber and the spin-atomic density as 



C a =Pa4>o/B- 



(24) 



This formula provides us the approach to measure 
whether the system is in the QSH phase. Firstly, we 
measured the spin- atomic density and denoted it as p2 , 
at p=0 in the absence of B. Then the optical lattice is ro- 
tated to generate the effective uniform magnetic field S, 
and the new density of the cold atoms p j ^ is measured 
again. If pi > p2 and p\ < p?, the system is in QSH 
insulator phase. However, if pi \=p+ i , the system is in 
the normal insulator phase. Since the density difference 
is actually quantized in units B/(f>o, the above method 
could be rather robust. 

It is clear from Eq. ([24]) that the total Chern number 
C=C^+Ci=0 and thus the direct detection of the QSH 
phase is a challenge in an electronic system. However the 
QSH phase can be directly verified in the atomic systems 
since the densities p^ and p± can be separately detected 
[5[ . To experimentally detect the spin-atomic density, we 
need first transfer the dressed state \xi) to a different 
hyperfine level \F) which is coupled to the excited state 
1 3) by a laser pulse (with a Rabi frequency Qf), as seen 
in Fig. 5(b). This pulse, together with the original laser 
beams Q p and Q c , make a Raman transition with an ef- 
fective Hamiltonian H R =(Q* F Q/A) \x±) (F|+H.c. (note 
that the state |xt) is still decoupled because of the phase 
relation between Q p and Q c ) [5]. Although the form of 
the state \xi) is spatially varying, the Rabi frequency Q 
(and thus also Q F Cl/A) is spatially constant. A complete 
Raman transition with a tt pulse will transfer all of the 



atoms being in the dressed state \xi) to hyperfine state 
\F). After this operation, the initial different dressed 
spin states are mapped to different hyperfine levels, and 
the populations in different atomic hyperfine levels can 
be separately imaged with the known experimental tech- 
niques. 



V. SUMMARY 

In summary, we have proposed a model which promises 
to host the transition from the QSH insulator phase to 
the normal insulator phase in the 2D Kagome optical 
lattice. The model includes two kind of periodic per- 
turbations, i.e., a nearest-neighbor intrinsic SO coupling 
and a trimerized Hamiltonian. The competition between 
them determines the system's phase. Then we demon- 
strate that the lattice SO coupling can be simulated by 
the laser-induced-gauge-field method and give the spe- 
cific laser setting and parameters to realize the intrinsic 
SO coupling. Furthermore, we have established the rela- 
tion between spin Chern number and spin-atomic density 
and then we can detect the spin Chern number through 
the standard density-profile technique used in atomic sys- 
tem. 
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